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ABSTRACT 



By combining //ewc/ie/-SPIRE data with archival Spitzer, H i , and CO maps, we investigate the spatial distribution of gas and dust in the two 

r\ ' famous grand-design spirals M99 and M 100 in the Virgo cluster. Thanks to the unique resolution and sensitivity of the Herschel-SFIRE pho- 

. , tometer, we are for the first time able to measure the distribution and extent of cool, submillimetre (submm)-emitting dust inside and beyond 

1-^ . the optical radius. We compare this with the radial variation in both the gas mass and the metallicity. Although we adopt a model-independent, 

Oh' phenomenological approach, our analysis provides important insights. We find the dust extending to at least the optical radius of the galaxy and 

showing breaks in its radial profiles at similar positions as the stellar distribution. The colour indices /350//500 and /250//350 decrease radially 

consistent with the temperature decreasing with radius. We also find evidence of an increasing gas to dust ratio with radius in the outer regions of 

both galaxies. 
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^O , 1- Introduction limited to the very nearby universe and large surveys of more 

C distant galaxies are unfeasible. 

t;;j- . The study of the gas and dust distribution in galaxies is essen- j . | _. 

in ■ tial to understanding their formation and evolution. The rate at The SPIRE instrument dGriffin et al.j |2010|) on-board 
^^ ; which gas is accreted and converted into stars regulates not only Herschel (jPilbratt et all |2010|) now bridges this gap observing in 
Q ■ the star formation history of galaxies but also their chemical evo- the range of 250/im-500//m. With the benefit of the stable condi- 
O lution. Dust is supposed to play a key role in this process. Dust tions of a space observatory, it is much more sensitive to the cold 
^H ! grains are the main coolant in star-forming galaxies, shielding component and provides excellent maps at high resolution, so is 
L; ■ the gas from the UV radiation and representing the site at which ideal for large surveys of many galaxies. With Herschel, we are 
. ^ ; H I is converted into H2 , and then collapses into stars (see e.g., "ow able to tackle the problem of the interplay between gas and 
Jv^ '. reviews by ibalzetti 2001' Draine 2003) dust, and combined with the number of recent hi gh resolution 
5h ■ To understand ihe dust, we ne^ed to' map the cold compo- surveys tracing the gas mass of local galaxies (e.g. |ChungetaH 
.^ ; nent, which does not dominate the energy but dominates the 2^9^ '^^""^t^l-' ^22^' ^^ '^^" ^tudy the distribution of gas, 
mass. However, using previously existing facilities, our knowl- "^^^^1^' ^^^ ^ust on a kpc scale for hundreds of galaxies, 
edge of the interplay between gas and dust and their radial distri- Here, we discuss first results for two famous grand-design 
bution have remained highly unce rtain, being based only onA< spirals M99 and M 100 (see Fig.[B in the Virgo cluster. We 
160iUm space obser vations (e.g. iMunoz-Mateos et al.L l2009at explore the distribution of the cool dust traced with Herschel 
iBendo et all 1201 Oal) . and challenging ground-based submm ob- by inspecting their radial profiles from mid-infrared to submm 
servations that were of optimal quality only at significantly wavelengths. We then attempt to correlate this with the ob- 
longer wavelength. For many galaxies, only integrated quanti- served gas and metallicity distributions and search for temper- 
ties have been derived because of the poor resolution of previ- ature variations. We use the Herschel-SPIRE maps taken dur- 
ous satelhtes. Resolved submm studies from the ground remain ing HerscheFs Science Demonstration Phase. T he two galax- 

ies are part of the Herschel Reference Survey (iBoselli et al.L 

* Herschel is an ESA space observatory with science instruments pro- |2010a|). This guaranteed time key project will provide maps for a 

vided by European-led Principal Investigator consortia and with impor- statistically-complete sam ple of 323 nearby galaxies i n all three 

tant participation from NASA SPIRE bands. In the RC3 dde Vaucouleurs et al.Lll991l) . the clas- 
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Fig. 1. The matched SDSS, MIPS 70 yum, SPIRE 250, 350, 500 yum, and H i -maps of M 99 f M/^per row) and M 100 f/owr row). The 
green contours are SPIRE 250 at levels about 0.03(a;5cr), 0.10, 0.34, 1.15, 3.9 Jy beam"'. The 350 and 500 //m are displayed at a 
high contrast to show the extent of the dust emission and the scan-artefact free smooth background dominated only by confusion 
noise. The white circle on the H i -maps marks the optical radius defined by D25. The blue circles on each map indicate the FWHM. 
Note, the SDSS jpg is scaled but not WCS matched. 



sification and optical radius R25 of M 99 and M 100 is given as 
SA(s)c, 2.69'(=« 12.9kpc), and SAB(s)bc, 3.71'(^ 17.8kpc), re- 
specti vely. We assume a distance for the Virgo Cluster of 16.5 
Mpc (lMeietal.Ll2007h . 



2. Data 



2.1. Herschel-SPIRE 



The SPIRE photometer dGriflin et al.L l2008l |20IO|) data were 
processed up to Level-1 (i.e., to the level where the pointed 
photometer time-lines were derived) with a customised pipeline 
script adapted from the official pipeline {POFS-pipeline.py, 
dated 27 Nov 2009) provided by the SPIRE Instrument Control 
Centre (ICCfl This Jython script was r un i n the Herschel 
Interactive Processing Environment (HIPE lOti l2OI0h coming 
with continuous integration build number 3.0.327, which is 
the current developer's branch of the data reduction software. 
However, in terms of the SPIRE scan-map pipeline up to Level- 
1, this is in principle identical to the Herschel Common Science 
System/Standard Product Generation v2. 1 , even down to the cal- 
ibration file^ associated with the individual pipeline modules. 
This version of the pipeline is used at ESA to produce the stan- 
dard products that will be available from the Herschel Science 
Archive once they become public. 

Currently, the Level-1 photometer time-lines still requires a 
residual baseline subtraction to be made. However, instead of 
subtracting the median of the time-line for each bolometer per 
scanleg (the default), we subtracted the median of the time-lines 
for each bolometer over the whole observation. This circumvents 
shadow artefacts caused in cases where the signal time-lines in 



' See 'The SPIRE Anal ogue Signal C hain and Photometer Detector 
Data Processing Pipeline' iGri ffin. 2009) for a more detailed description 
of the pipeline and a list of the individual modules. 

^ Aprt from the BsmPos file, for which we use an updated version 
that should improve the absolute astrometry. 



individual scanlegs are dominated by structured emission e.g. a 
large, extended galaxy or a strong cirrus component. 

This baseline-subtracted Level-1 data were then fed through 
an iterative de-striper, which minimises the difference between 
the signal in individual detector time-lines and the final map (see 
iBendo et al.. 2010b. for a longer description). At the end of this 
process, the signal time-lines were then mapped into a final im- 
age using the Naive Mapper available in HIPE. 

For M 99, the de-striping approach left some residual large- 
scale gradients. In this case, we resorted back to an initial base- 
line subtraction on a scan by scan basis. However, instead of a 
median, we used a robust linear fit with outlier rejection to the 
first and last fifty sample points, thus avoiding the galaxies in the 
centre of the time-lines. 

According to the IC C, the uncertain t y in th e flux calibration 
is of the order of 15% dSwinvard et al.Ll20IOl) and is currently 
based on a preliminary calibration. However, the ICC has re- 
leased some interim small correction factors to improve this cal- 
ibration. All flux values derived using the current standard cali- 
bration file for the flux conversion, are multiplied by 1.02, 1.05, 
and 0.94, for the 250yum, 350yum, and 500yum, respectivelyQ The 
full widths at half maximum (FWHM) of the SPIRE beams are 
18.1", 25.2", and 36.9", the pixel sizes are 6", 10", and 14" at 
250, 350, and 500 //m, respectively. 

For both, M99 and M 100, we observed a 12'xl2' field do- 
ing three repetitions of a cross-linked scan-map at nominal de- 
tector settings and nominal scan speed (30"/sec). The M 100 ob- 
servation was carried out twice. Both were treated independently 
here and used to verify the consistency of our results. 

2.2. MIPS, CO, and HI 

The 24, 70, and 160 ^m i mages were part of the SINGS sur- 
vey ( iKennicutt et al.L 2003") and were processed using the MIPS 
Data Analysis Tools (iGordon et al .. 2005) along with techniques 



See |http://herschel.esac.esa.inl/SDP-wkshops/presentations/IR/3-Griffin_SPIRE_SDP 
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Fig. 2. Left: Radial surface brightness profiles for M 99 (top) and M 100 (bottom) obtained from the smoothed and matched maps 
(see text). The profiles are from the bottom up: MIPS 24, 70, 160 ;um, SPIRE 250, 350, 500 yum (all in [log(Jy/arcsec2)]), H i (black) 
and H2 ^ {gray, both in [log(M0/pc^)]) adding 0.5, 0, 0, 1, 2, 3, -1, -1, respectively, to separate them in the plot and the radial axis 
is in units of 7?25. In all cases, we plot four different radial profiles per galaxy to show the influence of the fixed ellipse fitting (see 
Sect.O. Middle: SPIRE colours /350//500 and /250//350 with radius for M 99 (red) and M 100 (blue). Right: /350//500 and 
/70//160 versus the metallicity again for M 99 in red and M 100 in blue . Since the uncertainty in the SPIRE flux is still rather large, 
we only provide an indication of the error in the SPIRE colours with a black open circle in the lower left corner of the appropriate 
panels (see Sect. l2.3l l. The dash ed black line in the low er right panel is a fit to the integrated properties of a larger sample of galaxies 
with a variety of morphologies jBoselli et aUl20I0bl) . 



described bv iBendo et alJ (l20I0ah and Clements et al. (2010, 
in preparation). The FWHM of the MIPS beams are about 6", 
18", and 38", the pixel sizes are 1.5", 4.5", and 9" per pixel at 
24, 70, and 160 jum, respectively. The €0(7=1-0) maps, used as 
the tracer of the molecular hydrogen dominating the molecular 
mass, are taken from the No beyama CO Atlas of Nearby Spiral 
Galaxies (lKunoetalil2007h . The FWHM is 15" and the pixels 
are 1" per pixel. For the Hi, we used the zeroth moment maps 
from the VIVA survey (V LA Imaging of Virgo Spirals in Atomic 
Gas: IChungetaU 120091) . The FWHM and pixel sizes are ^ 37" 
and 5" per pixel for M 99 and ^ 30" and 10" per pixel for M 100. 

2.3. Analysis 

The SPIRE maps were first converted into Jy per pixel assum- 
ing a Gaussian beam (of the above quote d sizes). For t he Hi , 
we used the elliptical beam sizes given in IChung et al.l (l2009h . 
In all maps, we masked strong sources and artefacts. However, 
in the case of the SPIRE maps, being confusion limited, all faint 
sources are unmasked and part of the background. The residual 
background on each map was subsequently determined using 
IRAF'^ ELLIPSE task as described in Pohlen & Trujillo (2006). 
The background-subtracted maps of all wavelengths were there- 
after smoothed to the MIPS 160/im resolution of ^ 40" using 
custom convolution kernels derived as described in lBendo et alj 
(l2010ah and then matched to the SPIRE SOO/dm pixel size of 14". 
The SPIRE 250//m map was chosen to derive the final set 
of ellipse-fitting parameters (i.e., ellipticity, position angle, and 
centre). To ensure that our results are independent of the partic- 
ular ellipse geometry selected, we applied four different fixed 
ellipse fits to each map. For example, M 99 being a one-armed 



'' Image Reduction and Analysis Facility (IRAF) 
[http://iraf.noao. edu/| 



spiral, is slightly asymmetric (cf. Fig.[T]) so we selected one 
set of ellipse parameters derived in the outer parts and one in 
the inner parts, which each have slightly different centres. The 
final radial profiles, obtained using a combined mask on the 
smoothed and matched maps, are shown in Fig. [2] out to where 
we can trace signal on the map. The error bar in each measured 
point is a combination in quadrature of the uncertainty in the 
overall absolute calibration (currently for SPIRE the dominating 
source), the error in the ellipse intensity from the ellipse task, 
the uncertainty in the estimate of the background, and an 
additional uncertainty calculated by comparing the results from 
different versions of the pipeline. This last, very conservative 
uncertainty, is responsible for the currently rather large error 
bars in the measured flux ratios in Fig.|2] 



3. Results 

We detect dust emission traced by all three SPIRE bands out 
to at least the optical radius defined by R25 for both galax- 
ies (see Fig.[T] for ^ 3 - 5cr detections in the nominal maps, 
and Fig.|2] for the deeper, radially averaged profiles from the 
smoothed maps). Compared to H i , the dust can be found almost 
out to the Hi -edge of the regular disk for M 100. This is not 
entirely surprising, since M 1 00 is an intermediate H i -deficient 
(iHavnes & GiovanelliL IT981 ICavatte et al.L 1 1 994 galaxy (Hi- 
def= 0.35) and thus its outer Hi -disk has probably been al- 
ready stripped by the inte raction with the cluster environment 
(iBoselli & GavazziL l2006l) . A similar extension of the dust and 
H I -disk is observed for this range of Hi de ficiencies in other 
Virgo cluster galaxies (ICortese et al.L l2OI0l) . The situation is 
different for M99, which is not Hi -deficient (Hi-def= -0.1). 
Figure [T] clearly shows that the Hi emission is much more ex- 
tended than the submm at least to the north. Interestingly, this 
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extended H i halo (IChung et al.Ll2009l) however might be barely 
detected, but at the moment we cannot exclude that this is caused 
by residual background inhomogeneities coupled with a clus- 
ter of unresolved background sources. We can however exclude 
the presence of submm emission corresponding to the giant H i 
tail of M 99 (iHavnes et al.L l2007h to the southwest and we also 
find no measured flux associated with the extended l ow sur- 
face b rightness feature to the southwest of M 100 tChung et al.L 
I2OO9I) . 

The left column in Figure |2] shows the radial profiles 
obtained by the ellipse fitting to the smoothed maps. The 
MIPS and Hi profiles agr ee with the pub l ished ones by 
iMunoz-Mateos et al.1 (l2009bh and IChung et all (l2009l) . respec- 
tively. For M99, the MIPS and SPIRE profiles follow similar 
trends including a weak radial break in the profile, i.e., a change 
in the slope, at x 0.6 x R25, and a broken ex ponential is a more 
accurate fit than a single exponential (e.g. iPohlen & Truiillol 
I2OO6I for more background on b reaks). This break is also visi - 
ble in the optical profile shown bv lMunoz-Mateos et alJ (l2009bl) . 
The same is true for M 100, which also exhibits a more obvious 
break at around the same distance (it is ev en more striking in th e 
profile at native resolution presented bv ISauvage et al.L l2010h . 
This is the first time we see these breaks clearly in the dust distri- 
bution, while they are well-known at optical wavelengths. None 
of the so far presented hypotheses for the origin of these breaks 
have addressed this before (see e.g. the recent review by VlajiQ, 
I2OIOI for references) and it will be another pice of the puzzle to 
be explained by the various proposed models. The rising profile 
in the inner part of M 100 is related to the more prominent bulge, 
ba r, or inner-disk comp onent, which is discussed in more detail 
in lSauvageetal.l(l2010h . 

To investigate the variation in submm colours as a function 
of radius we plot in the middle column of Figure |2] the ratio of 
the SPIRE bands /350//500 to /250//350. These are colour 
temperature indices. The advantage of using these instead of 
a derived dust mass is that they are independent of the spe- 
cific, not yet well studied, model assumptions in this new wave- 
length range. They both decrease with radius, which suggests 
that the dust in the outermost regions is colder than in the cen- 
tre of the galaxies. This is naturally explained by an interstellar 
radiation field becoming less intense in the outski rts. Our pro- 
files are very similar to those of M 81 presented in lBendo et al.l 
( l2010bl) . who argue that the radial variation is driven by heat- 
ing from the evolved stars in the galaxy. The observed range of 
flux ratios along the galactic radii is the same as found for a 
sample of galaxies with a wide variety of m orphologies using 
integrated SPIRE fluxes (iBoselh et al.Ll2010bl) . Interestingly, the 
agreement between their integrated and our resolved analysis ex- 
tends beyond the colour profiles as shown in the right panels of 
Fig.|2] where we cou ple our colour gradien ts to the metallicity 
gradient published by SkiUman et al. ( 1996 ) renormalised to the 
[OIII]/[NII] base of iPettini & Pagell j20o4 . The trend we ob- 
serve radially for M 99 and M 100 matches the fit to the inte- 
grated properties of the lBoselli et al.l (1201 Obh sample very well. 
Both /350//500 and /70//160 (albeit only very weakly) de- 
crease with the radiaUy decreasing metallicity. This is again ex- 
pected since a l ower activity of star formation in the outer parts, 
as observed by IWilson et al.l (l2009t) . consequently entails lower 
metallicities. 

In Fig.|3]we finally show the ratio of the total gas mass (Hi 
plus H2 jj to 500yL(m flux ratio for the two galaxies. Since the 
500/im flux is a proxy of the dust mass, this provides a "model- 




We use a radially constant X-factor as given in lKuno et all ( 120071) 
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Fig. 3. Total gas mass (H i plus H2 ) to SPIRE 500//m flux ratio. 

independent" indication of the radial evolution of the dust-to- 
gas ratio. There is a clear trend visible with the gas-to-dust ra- 
tio increasing radially, which is consisten t with e arlier results 
(iMunoz-Mateos et all l2009al: iBendo et all l2010al) . but the ex- 
act shape needs to be revised once a proper SED dust modelling 
including the new SPIRE bands is available. 

In conclusion, we have found that the dust emission can be 
traced by the SPIRE bands at least out to the optical radius and 
beyond. The dust shows the same breaks in the radial profile as 
seen in the optical. The Hi is only slightly more extended but 
this needs to be regarded here in the context of the cluster envi- 
ronment. The SPIRE colour temperature indices decrease with 
radius following the measured trends in metallicity, and the ex- 
tent of the measured values along the galaxies' radii is consis- 
tent with the integrated properties of galaxies with a variety of 
morphologies. We have shown evidence of a radially rising gas- 
to-dust ratio. These results provide the first indication of the im- 
proved capabilities Herschel can offer for studying the resolved 
dust distribution in galaxies. 
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